A new model to describe the physics of (VO)2P207 
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In the past different models for the magnetic salt vanadyl pyrophosphate (VO)2P207 were discussed. 
Neither a spin ladder nor an alternating chain are capable to describe recently measured magnetic 
excitations. In this paper we propose a 2D model that fits better to experimental observations. 

PACS numbers: 75.I0.-b, 75.25.+Z, 75.40.Mg 



Low dimensional quantum spin systems have been a 
field of intense theoretical and experimental research over 
the last decades. Special interest was given to spin lad- 
der and chain materials. One compound that has been 
examined in this context is the insulating magnetic salt 
vanadyl pyrophosphate (VO)2P20y. Initially it was con- 
sidered as a prototypical realization of a two-leg antifer- 
romagnetic Heisenberg ladder [1]. However, susceptibil- 
ity data on polycrystalline and single crystalline material 
could be well fitted with both, ladder or alternating chain 
models [1—3], stressing the fact that susceptibility is not 
too sensitive to the particular model. Early inelastic neu- 
tron scattering measurements on polycrystalline samples 
indicated a spin gap of about 3.7 meV and supported a 
two-leg ladder model with the coupling constants esti- 
mated from susceptiblity data [4] . 

Recent neutron scattering experiments with powder 
samples [5] and with an array of single crystals [6] pro- 
vided detailed information on the low-energy excita- 
tion spectrum. Garrett et al. [6] observed a triplet 
branch with strongest (antifcrromagnetic) dispersion 
in b-direction, weak (ferromagnetic) dispersion in a- 
direction, and a spin gap of 3.1 meV. Most notably they 
found an additional second branch, separated from the 
first by an energy smaller than the gap. This was in- 
consistent with the picture, of (VO)2P207 being a spin 
ladder in a-direction, but also an alternating Heisen- 
berg chain in b-direction can not explain a second triplet 
branch over the whole Brillouin zone, as was shown re- 
cently [7,8]. 

In this work, starting with the alternating Heisenberg 
chain, we check whether coupling of (two) chains resolves 
this puzzling situation. As we do not succeed proceed- 
ing this way, we consider a new, truly two-dimensional 
model. We perform exact diagonalizations of finite sys- 
tems with up to 32 spins and periodic boundary condi- 
tions, supplemented by finite-size analysis if possible. 

The Hamiltonian of the alternating Heisenberg chain 
(AHC) reads as follows 



(1) 



where Si are spin-i operators and i denotes the sites in In- 
direction (see Fig. 1). For S > the spectrum has a gap; 
there is an one-magnon branch and a singlet branch, at 
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FIG. 1.: Schematic structure of (VO) 2 P 2 7 . The ex- 
change couplings are depicted for (i) the ladder model 



(J|| = J a ,J± = Jf, ), (ii) the alternating chain model 
(J* = J&(1 ± 5), and (iii) the new model (J*, J a , J* 
Throughout we measure energies in units of Jb- 



least around momentum ~, below a continuum of states. 

As an example Fig. 2 shows the low-lying excitations of 
a finite system of 32 sites for 5 = 0.2. The magnon branch 
is fitted to a sum of cosines ui™ = X)n=o a ™ cos (2«<z) and 
the shaded region corresponds to the continuum of two- 
magnon excitations resulting from this dispersion. Re- 
cently it was stressed [7,9] that there exists a second 
well-defined triplet below the two-magnon continuum 
near momentum but as Fig. 2 indicates, the second 
triplet occurs only very close to higher states, even for 
the relatively strong dimerization of 5 = 0.2. Therefore 
it was stated that an alternating chain will not explain 
the second triplet excitation observed in (VO)2P207 at 
all q-momenta. However, it is known [7] that including 
frustration, i.e. an antiferromagnetic next nearest neigh- 
bor interaction a between Sj and S i+2 j into the alter- 
nating chain model, yields a second well-defined triplet 
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FIG. 2.: Low-energy excitations of the AHC. The inset 
shows the difference A of the 1 st and 2 nd singlet and the 
2 nd triplet to the 1 st quintet at momentum 7r/2 versus 
inverse chain length. 

branch below the continuum in the whole Brillouin zone, 
provided a is sufficiently strong. 

Since an intra- chain frustration is not plausible in view 
on the structure of (VO)2P207, we will consider a per- 
pendicular coupling of two alternating chains instead: 

#cc = J b Yl 0- + <K- 1 ) i ) s ^- • s *+u 

3 = 1,2 

+ Ja ^ Si,l • Sj ; 2- (2) 

i 

Here j is numbering the chains. As was already suggested 
in [6] , such a coupling in a-direction should be ferromag- 
netic to explain the observed dispersion. For illustration, 
in Fig. 3 we plotted a few low-lying energies of a 2 x 12 
system with 5 = 0.2 and J a = —0.1. Similar results were 
obtained for a 3 x 8 system. 

Again we have a well-defined magnon branch. As a 
guide to the eye we shaded the region where one would 
expect a two-magnon continuum, approximated here by 
adding the gap at zero momentum to the magnon dis- 
persion. Close to the continuum edge there are several 
states: singlets, triplets, as well as quintets. To gain fur- 
ther insight one has to perform a finite-size analysis. For 
extrapolation to the infinite system we use the following 
formulas for the lowest singlets and triplets [7,10]: 

E S (L) = E s (oo) + (^ + C)e- L / A (3) 

E T (L)=E T (oc) + ^e- L / A (4) 

In Fig. 4 a few low-lying excitations at momenta 
{lb,Qa) = (0,0), (0, 7r) and (7r/2,0) are given subject 



FIG. 3.: Low-energy excitations of two alternating 
chains coupled ferromagnctically; the parameters are: 
chain-length = 12,8= 0.2, J a = -0.1 

to the ferromagnetic interchain coupling J a . At (0,0) 
and (0, 7r) the data is extrapolated to infinite system 
size, while at (7r/2,0) results for a 2 x 16 system are 
shown, because we have just four different system sizes 
(2 x 4, 8, 12, 16) at this momentum, making finite-size 
scaling questionable. Nevertheless, in the inset we tried 
to extrapolate the second triplet at (ir/2, 0) to the infinite 
system, using the ansatz of Eq. (4). The plot indicates 
that this triplet shows a weak nonmonotonic behaviour, 
in contrast to the single chain case. 

Obviously at momenta (0,0) and (0,7r) there arc no 
second triplets below the two-magnon continuum. Just 
at (qb,q a ) — ( 7r /2,0) a second triplet stays very close to 
the continuum edge, and the well-defined singlet excita- 
tion, known from the single alternating chain seems to 
disappear with increasing interchain coupling J a . 

From the above results we conclude that an interchain 
coupling of this simple type does not qualitatively change 
the structure of the low-energy excitations compared to 
the single alternating chain. Excitations are just shifted 
(as it seems linearly with J a in most cases), but no new 
features appear. This is, why we propose another model 
for (VO) 2 P 2 7 . 

We mentioned above that frustration in the alternating 
chain can lead to a well-defined triplet below the two- 
magnon continuum. Thus going to the second dimension 
we include an additional, frustrating coupling J x . Then 
our model Hamiltonian reads 

H x =J b Y / (l + 5{-l) i )S i j-S i+1J 

i,j 

+ Ja ^ Si,j ■ S,,j + 1 (5) 
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FIG. 4.: Energy of the lowest singlet (open circles) and 
triplet (filled squares) excitations for different interchain 
coupling J a and fixed dimcrization 5 = 0.2. 




FIG. 5.: Low-energy excitations of the 2D model; the 
system size is 4 x 8, 5 — 0.3, J a = 0.4 and J x = 0.425 



+ Jx / ,(S2i,j ' S2i+l,j+l + S2i+l,j ' S 2j j + i) 

i,j 

(cf. also Fig. 1, lower panel). As yet there is no data 
available about the strength of such a coupling, but as a 
first step it seems not unreasonable in view of the oxygen- 
mediated superexchange paths in (VO)2P207. We as- 
sume all exchange integrals to be antifcrromagnctic, but 
still the parameter space is very large. It appears that J x 
has to be bigger than J a to get a ferromagnetic magnon 
dispersion in a-direction, what is plausible. On the other 
hand, both couplings should not differ too much for a sec- 
ond triplet branch to exist in the whole momentum space, 
and should have a sufficient strength. The size of the gap 
to the first triplet branch is (still) mainly controlled by 
the dimcrization S. 

A good choice of parameters is 5 = 0.3, J a — 0.4 and 
J x = 0.425, for which we diagonalized systems of two, 
three and four chains with a total number of up to 32 
spins and periodic boundary conditions. The low-energy 
excitations of the 4x8 system are shown in Fig. 5. 

Beside two triplet branches we observe also a well- 
defined singlet, and there might even be a second singlet 
near momentum (qt,, q a ) = (tt/2, 0). As the difference be- 
tween J x and J a is small, the dispersion of the triplets is 
weak in a-direction, in accordance with experiments. We 
stress that the picture remains qualitatively unchanged 
going from the 3 x 8 to the 4x8 system, just the sec- 
ond triplet shifts downwards at momentum (0, 0) with 
increasing system size. Thus we believe that these fea- 
tures will survive in the infinite system. 

To provide some more information on the excitation 
spectrum, we calculated the dynamical spin structure 
factor and the integrated spectral weight 



S(q,co) = £ |HS*(q)|0)| 2 <5(£„ -Eo-u), (6) 

n 

N(q,w)= f dw'S(q,w'), (7) 
Jo 

where S z (q) = ^ • e' q Sfj . In the plot the integral 
is normalized to one; its real value N(q) = A(q, oo) is 
noted in each panel. 

It seems that we need a finite momentum component 
in b-direction for the first triplet to have some weight, 
while the second triplet appears only in a-direction. For 
comparison take the dashed and solid lines in the upper 
two panels of Fig. 6, corresponding to momenta (tt,x) 
and (0,x), respectively, that are equivalent in energy. 

To summarize, using exact diagonalization methods we 
have shown that a simple ferromagnetic coupling of al- 
ternating Hciscnbcrg chains does not provide two well- 
defined triplet branches as were observed in inelastic neu- 
tron scattering experiments on vanadyl pyrophosphate 
(VO)2?207. From our experience with frustrated al- 
ternating Heisenberg chains, we proposed an alternative 
model to describe the low-energy physics of (VO)2P207, 
introducing a frustrating interchain coupling. Due to the 
large parameter space and the computational effort for 
sufficiently extended 2D systems, we made no attempt 
to fix the parameters for (VO)2P207, but showed that 
the proposed model can describe the general feature of 
two triplet branches below a continuum of states. These 
triplets exhibit a ferromagnetic (antiferromagnetic) dis- 
persion in a- (b-) direction. Thus we believe that our 
model is a good starting point for further analysis. 

We thank E. Mullcr-Hartmann, J. Schlicmann, S. Sil 
and G. Wellein for stimulating discussion. G.B. ac- 
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FIG. 6.: Dynamical spin structure factor for the 2D 
model (4x8 system, S = 0.3, J a = 0.4, J x = 0.425); 
the lowest singlet and triplet excitations are classified. 
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